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ABSTRACT 

Experimental  and  theoretical  studies  have  continued  to  evaluate  the 
effects  of  surface  reactions  on  the  acoustic  response  function  of  composite 
solid  propellants.  Using  a  carboxy- terminated  polybutadiene-ammonium 
perchlorate  propellant,  the  effects  of  various  coatings  on  the  oxidizer 
crystals  were  investigated  using  a  T-burner  over  a  range  of  frequencies 
from  150  to  5,000  cps.  The  results  show  that  coatings  of  Kel-F  and  Viton 
reduce  the  acoustic  response,  the  reduction  increasing  with  increasing 
coating  level.  The  effect  results  from  the  coating  on  the  oxidizer  because 
incorporating  the  same  quantity  of  coating  material  in  the  binder  had  no 
affect  on  the  maximum  response. 

The  influence  of  both  chemical  structure  of  the  binder  and  curing  agent 
was  also  investigated  using  CTPB,  PBAA/AN,  and  UTREZ  binders,  and  MAPO  and 
epoxy  curing  systems.  The  data  show  PBAA/AN  propellants  have  approximately 
50%  of  the  acoustical  response  of  the  CTPB  propellants.  On  the  other  hand, 
the  cure  agent  was  responsible  for  greatest  shift  in  the  frequency  at  which 
the  maximum  admittance  occurs.  The  mechanism  responsible  for  this  frequency 
shift  currently  is  not  understood. 

Studies  were  made  of  the  effects  of  pressure  on  a  carboxy- terminated 
polybutadiene-ammonium  perchlorate  propellant.  The  results  indicate  that 
over  a  pressure  range  from  100  to  500  psia,  both  the  maximum  value  of  the 
acoustic  response  function  and  the  frequency  at  which  it  occurs  increase 
with  increasing  pressure. 

Theoretical  studies  were  conducted  using  a  simplified  combustion  model 
which  incorporates  the  effects  of  pressure -dependent  surface  reactions, 
pressure  dependent  gas  phase  reactions,  heat  transfer  and  surface  pyrolysis 
reactions.  Using  a  perturbation  approach,  the  predicted  acoustic  admittance- 
freqency  relation  is  characterized  by  three  independent  parameters.  Para¬ 
metric  studies  reveal  that  high  ratios  of  the  maximum  acoustic  response 
function  to  the  burning  rate  pressure  exponent  are  predicted  when  the  net 
heat  release  at  the  propellant  surface  (pressure-dependent  reactions  plus 
pyrolysis  reactions)  is  nearly  zero  or  is  exothermic.  The  stability  of 
the  combustion  process  in  self-excited  modes  has  also  been  considered 
theoretically  and  the  stability  bounds  determined  as  a  function  of  combus¬ 
tion  parameters.  The  perturbation  technique  was  expanded  to  examine  the 
nonlinear  behavior  of  the  combustion  model.  Preliminary  observations 
indicate  that  the  effect  of  oscillating  pressure  on  the  burning  rate 
changes  is  small  and  can  be  neglected. 
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a  «=  4  a  w  /r 

o 
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a^  *=  1  +  a 

A1  =  Mi(l+2Mr) 

1  2  2 
B  =  M  -M,  +  M 

r  i  r 

c  =  sonic  velocity  of  the  combustion  gases 
E1  a  Mr  +  4SMJ./X 


“  activation  energy  oi  exothermic  surface  reactions 
**  activation  energy  of  overall  pyrolysis  rate  at  surface 
■=  activation  energy  of  solid  phase  combustion  reactions 
*»  activation  energy  of  gas  phase  reactions 

D  frequency  of  the  acoustic  pressure  oscillations  at  the  point 

where  a  and  a,  measured 
g  d 

=  steady-state  heat  flux  from  gas  phase  combustion  zone  to  solid 
surface 


F  «  instantaneous  heat  flux  from  gas  phase  combustion  zone  to  solid 
surface 

Im(p/e)  =  imaginary  part  of  acoustic  response  function 
k  =  thermal  conductivity 
L  *=  length  of  combustion  chamber 
Mx  »  -1/2  (l+  Vl+T] 

m'*'  *»  order  of  gas  phase  reaction 
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=  steady-state  burning  rate  exponent 
«=  pressure  exponent  on  gas  phase  heat  flux 
a  instantaneous  pressure 

*  steady-state  pressure 

»  amplitude  of  the  acoustic  pressure  oscillation 
«  instantaneous  linear  burning  rate  of  propellant 
“  steady-state  linear  burning  rate  of  propellant 

«  g-rx 
♦ 

=  real  part  of  acoustic  response  function 
**  La  Placian  variable 
■  1^0 

**  time 

•  temperature 

«  steady  state  temperature 

■  propellant  flame  temperature 

“  initial  temperature  of  the  propellant 
=»  distance  into  the  solid  from  the  surface 

»  Arrhenius  pre- exponential  factor  for  exothermic  surface  reations 

■  Arrhenius  pre -exponential  factor  for  overall  pyrolysis  rate  at 
surface 

=  Arrhenius  pre -exponential  factor  solid  phase  combustion  reactions 
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Greek  Symbols 

propellant  thermal  diffusivity 

dP^/dt  during  the  decay  of  the  acoustic  pressure  oscillations 

dP^/dt  during  the  growth  of  the  acoustic  pressure  oscillations 

AH  /C  T°(0) 
v  p 

Pq  Z  exp(-E1^T°(0))/fop  CpT  (0) 

effective  heat  of  gasification  of  propellant  at  the  surface 
perturbation  variable 

amplitude  of  the  first  harmonic  pressure  oscillations 
amplitude  of  second  order  pressure  oscillations 

instantaneous  value  of  ^ 

r  x 
o 

a 

F  /r  oC  T°(0) 
o  o  p 

component  of  temperature  oscillating  with  first  harmonic 

component  of  temperature  oscillating  with  second  harmonic 

steady  state  state  shift  in  temperature  resulting  from 
second 

amplitude  of  second  harmonic  temperature  oscillations 
resulting  from  second 

VE2 

ka  w/r 

o 

4a  w/Z22 

first  harmonic  amplitude  of  H- 

first  harmonic  oscillating  component  of  burning  rate 
(E2  0(O,t)/RT°-(O) 

0o2(O)/T°(O)^ 
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p  -  propellant  density 
T  =  dimensionless  time 

<f>  =  (0) 

=  RT°  (0)/E2 
w  =  2tt  f 
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> 1  =  complex  frequency  = 

v  =  0.5  tan"1  (X.) 
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1.0  INTRODUCTION 

Unstable  combustion  of  solid  rocket  propellants  is  of  interest  for  two 
reasons:  first,  it  is  one  of  several  methods  available  for  investigating 
the  structure  of  the  solid  propellant  combustion  and,  second,  unstable 
combustion  has  presented  serious  problems  in  the  development  of  some  opera¬ 
tional  solid  propellant  rocket  motors.  For  these  reasons,  UTC  has  been 
investigating  certain  aspects  of  this  problem,  namely,  pressure-coupled 
combustion  instability,  under  Contract  No.  DA-04-200-AMC-968(X) .  The 
principal  objective  of  these  studies  is  to  examine  the  relation  between 
combustion  instability  as  observed  by  pressure  oscillations  in  the  gas 
phase  and  exothermic  reactions  occurring  on  and  within  the  solid  phase 
surface. 

Experimental  studies  have  been  direcced  toward  investigating  the 
effects  of  coatings  on  the  solid  oxidizer  crystals;  effects  of  the  chemical 
and  physical  nature  of  the  propellant  binder,  catalysts,  and  cross-linking 
agents  on  the  acoustic  admittance;  nonlinear  response  properties;  and 
propellant  burning  rate.  Concurrent  theoretical  studies  have  been  directed 
toward  predicting  the  effects  of  these  reactions  on  the  acoustic  response 
function.  These  studies  have  included  the  study  of  the  influence  of 
pressure-dependent  heat  transfer  from  the  gas  phase  combustion  zone,  the 
effects  of  various  combustion  parameters,  and  the  investigation  of  sub¬ 
surface  reaction  contributions.  This  report  presents  the  results  which  have 
been  obtained  during  the  second  year  of  study. 
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2.0  TECHNICAL  DISCUSSIONS 

Recent  studies  on  the  structure  of  the  combustion  zone  of  composite 
solid  propellants  have  presented  indications  that  there  are  exothermic  proc¬ 
esses  which  occur  on  and  within  the  solid  phase. ^  through  6 )*  The  results 
of  these  studies  indicate  that  the  combustion  process  is  controlled  by  two 
interdependent  exothermic  reaction  zones  near  and  on  the  surface  of  the  pro¬ 
pellant.  One  zone  is  in  the  gas  phase  at  a  finite  distance  from  the  solid 
propellant  surface  and  is  characterized  by  i-'ferdif fusion  of  gasified  oxi¬ 
dizer  and  fuel  species  and  combustion  of  particles  of  ejected  matter  from 
the  surface.  The  second  reaction  zone  occurs  on  and  within  the  solid  pro¬ 
pellant  surface.  The  primary  heat  release  in  this  zone  can  occur  as  a 
result  of  exothermic  decomposition  processes  and  chemical  reactions  between 
the  initial  decomposition  products  of  the  solid  oxidizer  and  the  adjacent 
fuel  surface.  Transient  and  steady-state  combustion  studies  indicate  that 
much  of  the  pressure-dependent  combustion  process  can -be  associated  with 
these  interfacial  reactions. 

The  exothermic  surface  reactions  at  and  within  the  surface  release 
sufficient  heat  to  expel  partially  combusted  products,  pyrolysis  products, 
and  fuel  and  oxidizer  fragments  into  the  gas  phase  zone  above  the  surface, 
where  they  intermix  and  burn  completely.  The  maximum  flame  temperature  is 
reached  in  the  luminous  zone,  where  the  largest  portion  of  the  heat  is 
released.  However,  because  of  the  relatively  large  mass  flow  perpendicular 
to  the  surface,  only  a  small  amount  of  heat  released  in  the  luminous  flame 
zone  can  normally  reach  the  surface  to  supplement  the  heat  generated  by  the 
interfacial  reactions. 

Incorporation  of  the  exothermic  chemical  processes  on  and  within  the 
solid  phase  represents  an  important  addition  to  the  analysis  of  propellant 
combustion  phenomena.  Previous  theoretical  treatments  of  steady-state 
combustion*-' *0'  as  well  as  combustion  instability  have  con¬ 

sidered  the  exothermic  combustion  reactions  to  occur  only  in  the  gas  phase. 
However,  Friedly (13, 14,15)  and  Capener(^)  have  shown  theoretically  that 
the  principal  time  lags  in  the  combustion  zone  are  associated  with  energy 
transport  within  the  solid  phase  for  oscillations  below  10,000  cps.  By 
comparison,  the  time  constants  associated  with  the  reaction  and  transport 
processes  in  the  gas  phase  are  small  and  these  processes  can  be  considered 
to  be  in  equilibrium  at  any  instant.  Because  the  reactions  on  and  below 
the  surface  are  exponentially  dependent  on  the  solid  phase  temperature,  the 
energy  release  by  these  reactions  should  be  an  important  consideration  in 
the  analysis  of  transient  combustion  phenomena.  Furthermore,  the  reported 
data  indicate  that  these  reactions  on  and  within  the  solid  phase  are  first 
order  in  local  pressure.*-^'  Thus,  the  acoustic  pressure  oscillations  will 
produce  oscillations  in  reaction  rates  through  both  the  pressure  and 
temperature  sensitivities  of  the  reaction  rates  which  can  have  a  signifi¬ 
cant  effect  on  the  combustion  stability. 

*  Parenthetical  superscript  numbers  denote  references  appearing  on  page  42. 
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2.1  EXPERIMENTAL  STUDIES 

The  experimental  studies  conducted  as  part  of  the  research  program 
were  directed  toward  investigating  the  effects  of  propellant  combustion 
variables  on  the  response  of  the  combustion  process  to  acoustic  pressure 
oscillations.  The  principal  propellant  property  derived  from  these  studies 
was  the  real  part  of  the  acoustic  response  function,  Re  (fi/e),  which  relates 
the  amplitude  of  the  burning  rate  oscillations  to  the  amplitude  of  the 
incident  acoustic  pressure  oscillations  by  the  relation 


Re 


(«)■ 


(rTfo) 

(?/Po) 


(1) 


2.1.1  Experimental  Apparatus 

Measurements  of  the  acoustic  response  function  were  obtained  in  a  T- 
burner  apparatus.  This  type  of  equipment,  which  has  become  almost  standard 
for  acoustic  measurements,  consisted  of  a  cylindrical  combustion  chamber 
with  end -burning  charges  of  propellant  in  either  or  both  ends.  The  appara¬ 
tus  used  in  these  studies  was  based  on  the  design  concepts  of  Horton^') 
and  Stittmater and  has  an  internal  diameter  of  1.5  in.  over  the  entire 
length.  Propellant  samples  are  located  in  end  caps  which  also  contain  pro¬ 
vision  for  the  igniter  leads  and  pressure  transducers.  Electrically  heated 
nichrome  ignition  wires  are  soldered  to  two  electrodes  which,  in  turn,  are 
located  in  a  ceramic  tube  inside  the  cap  and  pressure  sealed  with  a  connax 
fitting.  The  nichrome  wire  is  passed  through  the  propellant  sample  and  is 
secured  to  the  surface  with  small  staples.  Powdered  propellant  is  then 
placed  over  the  wires.  Combustion  products  and  pressurized  gas  are 
exhausted  into  two  blowdown  tanks  through  four  orifices  equally  spaced 
around  the  circumference.  Details  of  the  apparatus,  experimental  procedure, 
and  the  electronic  data  acquisition  system  have  been  reported  previously. 


The  relationship  between  the  observed  acoustic  pressure  oscillations 
and  the  acoustic  response  function  of  the  propellant  has  received  consider¬ 
able  attention.^' >18,20)  Although  several  methods  have  been  developed  to 
derive  the  acoustic  response  function  from  acoustic  pressure  measurements, 
the  only  method  used  in  this  study  was  based  on  growth  and  decay  constant 
determine .ions .  Hence,  the  expression  used  to  reduce  the  oscillating 
pressure  measurement  to  the  acoustic  response  function  is 


Re 


4c  p  r 
s  o 


a 

_£  .. 


01, 


(2) 


During  the  experimental  program,  special  care  was  required  when  tests 
were  conducted  at  low  frequencies.  Under  these  conditions,  the  combustion 
chamber  becomes  long  and  the  quantity  of  room  temperature  gas  required  to 
pressurize  the  chamber  thereby  becomes  significant.  After  ignition,  a 
significant  period  of  time  is  required  to  achieve  thermal  equilibrium  in 
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the  combustor.  However,  the  frequency  of  the  pressure  oscillations  depends 
on  the  temperature  through  the  sonic  velocity,  i.e,,  c~  >/T7  and  the 
expression 


f  =  c/2L  (3) 

8 

One  point  worthy  of  note  is  that  the  sonic  velocity  which  appears  in 
equation  2  results  from  considerations  at  the  combustion  zone  and  therefore 
should  be  evaluated  from  the  experimental  flame  temperature.  The  approach 
used  in  this  study  is  based  on  the  assumption  that  at  high  frequencies 
(short  combustion  chamber),  the  experimental  flame  temperature  is  established 
throughout  the  combustion  chamber  very  rapidly.  Thus,  c  was  evaluated  from 
equation  3  by  extrapolating  the  parameter  2Lfg  to  high  frequencies 

2.1.2  Propellant  Formulations 


There  are  several  potential  methods  for  controlling  the  contributions 
of  surface  reactions  to  the  acoustical  response.  One  method  involves 
coating  the  oxidizer  crystals  with  coating  materials  in  an  attempt  to  alter 
the  reactivity  of  the  crystal-binder  interface.  This  approach,  which  was 
examined  during  the  first  year  of  effort,  has  been  extended  during  the 
second  year.  The  details  of  the  coating  and  propellant  preparation  have 
been  presented  by  Brown,  et  al.(^) 

In  addition  to  the  effects  of  coatings,  changes  in  the  chemical 
structure  of  the  binder  and  in  the  agents  used  to  crosslink  the  binder  were 
suspected  of  having  a  significant  influence  on  the  acoustic  response  func¬ 
tion.  Accordingly,  propellants  were  prepared  which  permitted  investigation 
of  these  variables. 

The  particular  propellant  formulations  tested  and  the  resultant 
burning  rate  data  are  shown  in  table  I.  Burning  rate  measurements  were 
made  in  a  standard  strand  burner  at  three  pressure  levels:  200,  500,  and 
1,000  psig.  The  data  for  each  system  were  checked  at  the  200-psig  level 
with  the  total  burning  time  measurements  in  the  T-burner,  i.e.,  sample 
size  with  1/4- in. -thick  samples.  The  difference  between  the  two  values  was 
less  than  10%  in  all  cases.  Because  of  the  experimental  scatter  in  the 
burning  rate  measurements,  these  data  are  valid  to  within  ±10%  of  the 
reported  values.  Also,  the  reported  burning  rate  pressure  exponents  are 
abnormally  high,  especially  considering  the  large  oxidizer  particle  size 
which  was  used.  However,  the  data  at  higher  pressures  show  the  exponents 
for  many  of  these  propellants  are  not  constant  over  the  entire  range. 

2.1.3  Experimental  Results 

During  the  second  year  of  study,  the  acoustic  response  function  of 
propellants  10  through  15  were  determined  as  a  function  of  frequency  in  the 
T-burner.  These  data  are  presented  in  figures  1  through  7  using  the 
normalized  frequency  8  x  fg/r^  as  the  parameter.  This  form  was  used 
based  on  the  theoretical  studies  of  Hart^H'  and  Friedly(13-15)  and  the 
theoretical  studies  to  be  discussed  in  the  sections  that  follow.  These 


4 


UTC  2136-TSR2 


theoretical  studies  suggest  that  the  proper  form  for  the  dimensionless 
frequency  should  include  the  propellant  thermal  diffusivity.  For  the  pro¬ 
pellant  listed,  this  propellant  property  should  be  relatively  constant 
(approximately  3x10"^  in.^/sec)  and  therefore  has  been  omitted.  In 
addition,  the  definition  of  the  acoustic  response  function  is  such  that  it 
should  approach  the  normal  burning  rate  pressure  exponent  as  the  frequency 
approaches  zero;  therefore,  the  data  are  reported  showing  the  extrapolation 
at  low  frequencies  to  the  burning  rate  pressure  exponent  obtained  from  the 
burning  rate  measurements. 

For  completeness,  all  the  experimental  data  obtained  for  all  the  pro¬ 
pellants  shown  in  table  I  are  presented  in  appendix  A.  The  values  of  the 
growth  frequency  (fg),  growth  rate  (#g),  decay  frequency  (f,),  decay  rate 
--(o'd)  and  Re(p/e)  are  reported  with  the  graphs  of  the  data  snowing  the 
acoustic  response  as  a  function  of  the  measured  growth  frequency. 

It  is  important  to  note  that  the  data  are  subject  to  significant  un¬ 
certainty.  The  steady-state  burning  rates,  as  previously  indicated,  may 
have  errors  as  high  as  ±10%.  These  values  enter  into  the  calculation  of 
the  response  functions,  i.e.,  u  oc  l/fQ,  and  the  frequency  functions,  8ir  fg/ 
r0^.  In  addition  to  these  uncertainties,  the  measurements  of  the  growth 
and  decay  rates  are  subject  to  an  uncertainty  of  ±10%;  therefore,  the 
absolute  value  of  the  response  function  u/e  is  valid  to  within  ±20%.  The 
experimental  results  must  be  interpreted  accordingly. 

The  data  obtained  using  the  propellants  containing  3%  Kel-F  and  1.5% 
Viton  A- coated  AP  are  shown  in  figure  1,  together  with  previously  reported 
data  on  the  1.5%  Kel-F-coated  AP  and  the  1.2%  Kel-F  in  the  binder.  The 
effect  of  these  coatings  on  the  maximum  admittance  is  significant.  Coating 
the  oxidizer  with  Kel-F  800  reduced  the  maximum  admittance  by  50%  compared 
to  the  uncoated  oxidizer.  Increasing  the  level  of  coating  further  reduced 
the  response  function,  although  the  reduction  does  not  increase  linearly 
with  coating  level.  The  effect  of  the  Kel-F  coating  on  the  magnitude  of 
u/e  is  consistent  with  the  concept  of  interfacial  reactions  because  Kel-F 
is  more  resistant  to  oxidation  than  the  CTPB  binder.  However,  the  data 
reported  previously(^)  on  the  effects  of  other  coatings  (Hypalon  30  and 
ethylcellulose)  showed  no  effect  of  the  coating  on  the  maximum  acoustic 
response  function.  Apparently,  only  the  Kel-F  and  Viton  reduced  the 
reactivity  of  the  exothermic  reactions  in  the  solid  phase  sufficiently  to 
have  a  significant  effect  on  the  response  properties. 

Because  both  these  materials  contain  halogen,  one  is  tempted  to  con¬ 
clude  that  the  presence  of  halogen  is  responsible  for  the  changed  behavior. 
However,  the  data  indicate  that  incorporating  the  Kel-F  into  the  binder 
appears  to  have  had  very  little  effect  on  the  response  function.  The 
observed  effect  is  approximately  a  10%  increase,  but  this  effect  is  less 
than  the  experimental  error.  Hence,  the  location  of  the  additive  is  also 
important  and  it  is  postulated  that  the  contribution  of  interfacial  reac¬ 
tions  to  the  overall  combustion  process  is  in  part  reduced  by  the  coating, 
which  in  turn  reduces  the  admittance. 
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Also,  it  is  interesting  to  note  that  the  burning  rate  and  the  burning 
rate  pressure  exponent  for  the  standard  propellant,  as  well  as  the  Kel-F- 
modified  and  Kel-F-coated  AP  essentially  are  identical.  When  compared  with 
the  acoustic  response  data,  these  results  would  suggest  that  the  postulated 
effect  of  interfacial  reactions  is  more  significant  under  transient  combus¬ 
tion  conditions  than  under  steady- state  conditions. 

Another  aspect  of  the  experimental  studies  was  concerned  with  the 
large  difference  between  the  results  obtained  by  Horton^^'  and  the  pro¬ 
pellants  tested  in  this  study.  The  maximum  acoustic  response  of  about  6 
was  found  for  the  CTPB  propellants,  where  Horton  observed  the  maximum  of 
about  3  for  the  PBAA  propellant.  It  appears  that  the  binder  has  a  signifi¬ 
cant  effect  on  the  acoustic  response  function.  To  study  this  effect  in 
greater  detail,  propellants  were  prepared  in  which  the  CTPB  binder  was  re¬ 
placed  by  PBAN  cured  with  MAPO  in  one  case  and  with  epoxides  in  another. 

The  results  are  shown  in  figure  2  for. both  the  MAPO  cross-linked  binder 
(UTX-9162)  and  the  epoxy  cross-linked  propellant  (UTX-9168).  The  epoxy- 
cured  propellant  had  a  burning  rate  at  200  psig  of  r  =  0.145  in. /sec  and 
the  MAPO-cured  propellant  of  r  =  0.127  in. /sec.  For  comparison,  the  data 
for  the  MAPO-cured  CTPB  propellant  are  also  shown  in  figure  2  and  indicate 
that  a  change  of  binder  reduced  the  acoustic  response  significantly. 

Because  the  magnitude  of  the  response  function  measured  by  Horton  on  PBAA 
propellants  and  these  data  are  nearly  equal,  these  results  suggest  that  the 
data  appear  to  be  consistent  between  the  two  T-burners.  It  would  appear 
that  the  choice  of  tne  CTPB  binder  was  particularly  fortunate  for  studying 
the  additive  effects.  It  also  appears  that  CTPB  propellants  have  somewhat 
higher  acoustic  response  functions,  at  least  for  nonaluminized  AP  formula¬ 
tions  loaded  to  75%  to  78%. 

The  effect  of  changing  the  MAPO  curing  agent  used  in  the  previous 
studies (19)  to  an  epoxy  system  has  been  also  studied  using  the  standard 
CTPB  propellant  formulation  with  noncoated  oxidizer.  The  response  functions 
for  UTX-9173  are  shown  in  figure  3,  which  was  an  epoxy-cured  CTPB  propellant 
(78%  AP),  along  with  the  basic  formulation  UTX-8501  (78%  AP  in  CTPB  binder- 
MAPO  cure).  The  epoxy-cured  propellant  had  a  burning  rate  at  200  psig  of 
r  =  0.138  in. /sec  and  the  MAPO-cured  propellant  was  r  =  0.113  in. /sec.  The 
results  indicate  that  the  position  of  maximum  acoustic  response  occurs  at  a 
lower  frequency  for  the  epoxy-cured  propellants.  However,  the  maximum  value 
of  the  acoustic  response  function  is  not  significantly  altered  by  the  change 
in  the  cure  system. 

Two  observations  can  be  made  from  these  results.  First,  the  acoustic 
admittance  for  the  PBAN  propellants  are  significantly  less  than  those  of 
the  equivalent  CTPB  propellant.  Secondly,  the  admittance  of  the  epoxy- 
cured  propellant  is  less  than  the  admittance  for  the  MAPO-cured  PBAN 
propellant.  This  latter  effect  was  not  observed  with  the  CTPB  propellant. 
Thus,  there  appears  to  be  an  interaction  between  the  cure  agent  and  the 
binder  such  that  cure  agent  effects  depend  on  the  binder.  The  mechanisms 
responsible  for  this  interaction  are  not  immediately  apparent  and  further 
study  is  required  to  explain  these  observations. 
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Figure  2.  Effect  of  PBAN  Binders 
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The  effect  of  ambient  pressure  on  the  acoustic  admittance  function  was 
also  studied  using  the  basic  78%  AP-CTPB  binder  propellant.  Pressures  of 
100,  200,  and  500  psig  were  investigated  and  the  results  are  shown  in  figure 
4.  The  dimensional  frequency  8  jt  f  /r  2  is  based  on  the  burning  rate  of 
the  propellant  at  the  specified  prelsure.  The  data  indicate  that,  over  the 
range  of  pressures  studied,  the  maximum  value  of  Re(p/e)  increased  as  the 
pressure  increased  (i.e.,  compare  200  and  500  psig).  This  latter  effect  is 
caused  by  r2  factor  in  the  dimensionless  frequency  (8  rf  f_/r2). 

a 

To  compare  these  results  with  pressure  data  obtained  on  a  76%  AP  (80u) : 
24%  PBAA  propellant  by  NOTS,(22)  the  reSp0nse  function  at  the  various 
pressure  levels  is  shown  as  a  function  of  the  measured  growth  frequency  in 
figure  5.  The  results  are  consistent  with  the  observations  by  NOTS  in  that 
the  greatest  shift  in  the  frequency  at  which  the  maximum  appears  in  the 
response  function  occurs  between  100  and  200  psig  and  only  a  small  change 
occurs  between  200  and  500  psig 

Data  were  obtained  on  NOTS  propellant  A-13  in  connection  with  the  round 
robin  for  the  Standardization  Committee  on  Combustion  Instability  Measure¬ 
ments  and  was  presented  at  the  committee  meeting  on  6  June  1967  at  Anaheim, 
California.  The  data  were  obtained  over  a  frequency  range  from  250  to 
4,000  cps  at  a  chamber  pressure  of  200  psig  and  are  tabulated  in  appendix 
A.  A  comparison  of  these  results  with  the  earlier  work  of  Horton'21'  is 
shown  in  figure  6.  The  results  agree  favorably  over  this  frequency  range. 

In  the  frequency  range  from  1,000  to  4,000  cps,  the  present  results  were 
also  consistent  with  the  majority  of  investigators  on  the  committee. 

Further  experiments  were  conducted  with  a  84%  AP-UTREZ  (polyisobuty¬ 
lene)  binder.  (This  propellant  was  prepared  for  the  photographic  studies 
discussed  in  section  2.1.4.)  The  acoustic  response  data  for  the  UTREZ 
propellant  are  shown  in  figure  7  with  the  data  from  the  standard  847»  AP- 
CTPB  binder;  both  propellants  were  cured  with  the  same  curing  agent.  The 
burning  rate  of  the  UTREZ  propellant  was  0.125  in. /sec,  and  the  CTPB 
propellant  had  a  value  of  0.113  in. /sec  at  200  psia.  However,  both  the 
oxidizer  loading  as  well  as  the  binder  type  were  varied.  Current  studies 
include  measurement  on  a  847.  AP- loaded  UTREZ  propellant  and  the  signifi¬ 
cance  of  the  available  data  cannot  be  determined  until  further  data  can  be 
obtained. 

2.1.4  Photographic  Studies 

The  theoretical  model  for  the  combustion  process  depends  on  a  know¬ 
ledge  of  the  chemical  and  fluid  mechanical  processes  which  influence  the 
acoustic  admittance.  The  previously  discussed  T-burner  experiments  have 
been  useful  in  determining  the  effects  of  propellant  changes  on  the 
measured  pressure- time  variations  (i.e.,  acoustic  admittance),  but  it  is 
difficult  to  resolve  the  details  of  the  combustion  behavior  from  these 
types  of  indirect  observations.  Therefore,  it  is  desirable  to  examine  the 
propellant  when  it  is  burning  under  oscillatory  conditions  in  the  T-burner. 

A  sample  holder  was  constructed  with  viewing  windows  as  shown 
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Figure  4.  Effect  of  Pressure  on  Acoustic  Response 


Figure  5.  Effect  of  Pressure  on  Acoustic  Response 
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schematically  in  figure  8.  The  1-1/2-in. -diameter  propellant  disk  is 
located  in  the  center  and  a  purge  system  is  placed  between  the  sample  and 
the  windows.  Tests  were  conducted  using  both  Vycor  and  Plexiglas  windows 
with  and  without  the  purge  system.  The  windows  were  sealed  to  their 
mounting  with  gaskets.  A  photograph  of  the  final  assembly  is  shown  in 
figure  9. 

The  camera  used  in  these  studies  is  a  Hycan  16-mm  rotary  prism  camera 
with  a  framing  rate  capacity  to  4,000  frames/ sec,  which  will  provide 
several  frames  during  a  complete  cycle  at  the  lower  frequencies.  The 
initial  photographs  were  taken  at  a  framing  rate  of  500  frames/sec.  A 
quartz- iodine  lamp  was  used  for  backlighting  the  viewing  section.  The 
complete  camera  arrangement  is  shown  in  figure  10. 

The  greatest  problem  encountered  in  the  photographic  studies  is  the 
rapid  buildup  of  carbonaceous  material  on  the  windows.  To  eliminate  this 
problem,  several  changes  were  made  in  the  experiments.  A  "clean  burning" 
propellant  consisting  of  84%  AP  was  selected.  The  burning  rate  at  200  psig 
is  0.130  in. /sec.  The  windows  were  also  recessed  from  the  tube  walls  to 
allow  the  nitrogen  window  purge  to  operate  efficiently.  The  amount  of 
ground  propellant  placed  ever  the  ignition  wires  was  reduced  to  minimize 
the  effects  of  the  additional  mass  ejected  during  the  ignition  phase  and 
subsequent  clouding  of  the  windows. 

These  changes  significantly  improved  the  initial  photographic  results 
and  made  it  possible  to  observe  the  propellant  surface  during  a  large  por¬ 
tion  of  the  run.  During  the  early  stages,  a  film  of  opaque  material  is 
temporarily  deposited  on  the  windows,  thereby  preventing  observation  of 
the  initial  stages  of  the  instability  process.  Methods  of  solving  this 
problem  are  currently  being  considered. 

A  typical  film  sequence  is  shown  in  figure  11  for  the  84%  AP 
propellant.  The  sides  of  the  test  sample  were  inhibited  before  firing  and 
a  nitrogen  purge  was  used  in  this  test.  The  photographs  were  taken  from 
the  16mm  color  film,  and  the  original  film  displays  greater  detail  than 
appears  in  figure  11.  Random  frames  were  selected  to  display  various 
features  of  the  combustion  process.  The  first  frame  shows  the  test  section 
before  firing  and  the  next  frame  displays  tha  ignition  process.  Ignition 
is  obtained  by  electrically  heating  the  nichrorae  wire  secured  to  the  pro¬ 
pellant  surface  and  this  appears  as  a  blue  flash  (i.e.,  bright  flash  near 
surface  in  frame  2)  when  the  electrical  contact  is  broken.  Once  ignition 
has  occurred,  mass  is  evolved  from  the  surface  as  displayed  by  the 
luminous  streaks  near  the  surface  in  frame  3.  At  a  later  time,  the  pro¬ 
pellant  started  to  burn  around  the  inhibitor,  as  shown  in  frames  4  and  5; 
this  particular  problem  is  presently  under  investigation. 
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Figure  8.  Schematic  of  Viewing  Section 
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Figure  9.  Photograph  of  Viewing  Section 


Figure  11.  Photograph  of  Combustion  Process  in  T- Burner  Viewing  Secti 
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2.2  THEORETICAL  STUDIES 

The  theoretical  study  of  the  approximate  combustion  model  reported  in 
the  first  technical  summary  report  (UTC  2136-TSR1)  has  been  continued  in 
order  to  improve  the  basis  for  interpretation  and  correlation  of  the  exper¬ 
imental  data.  These  studies  were  motivated  by  two  considerations.  First, 
previously  reported  theories  have  not  accounted  for  the  magnitude  of  the 
response  functions  found  experimentally  in  this  program  and  those  reported 
by  Horton.^1'  Second,  earlier  studies  have  not  considered  the  effects  of 
exothermic  surface  reactions  of  the  type  discussed  in  preceding  paragraphs. 
Without  a  suitable  theoretical  framework,  it  has  not  been  possible  to 
determine  any  contribution  of  exothermic  surface  reactions  on  the 
acoustic  response  function.  Therefore,  development  of  a  model  which 
included  surface  reaction  effects  was  initiated. 

The  approach  taken  in  these  studies  has  been  to  formulate  a  simplified 
analysis  to  facilitate  the  study  of  the  overall  linear  behavior  of  the  model. 

In  addition,  these  studies  will  permit  estimation  of  the  important  non- 
linearities,  i.e.,  burning  rate  charges,  sinusoidal  wave  forms,  and  variation 
of  admittance  with  pressure  amplitude.  The  analysis  is  intended  to  be  expanded 
at  a  later  date,  as  dictated  by  the  experimental  data  and  the  results  of  the 
preliminary  theoretical  analysis,  to  study  the  more  important  processes  in 
detail. 

During  this  past  year  of  study,  the  basis  for  the  model  has  been  refined, 
the  linear  solutions  have  been  studied  in  greater  detail,  and  the  nonlinear 
behavior  has  been  studied  further.  The  following  paragraphs  present  new 
developments  in  these  areas. 


2.2.1  Analytical  Combustion  Model 

From  the  description  of  the  combustion  zone,  the  transient  behavior  of 
the  solid  phase  can  be  modeled,  at  least  to  a  first  approximation,  by  the 
one-dimensional  form  of  the  transient  heat  conduction  equation  expanded  to 
include  the  effects  of  exothermic  reactions.  This  equation  is 


»  a  +  r  +  Z  exp  (-E  /RT)  (4) 

St  0  %  Sx  3  3 

At  the  surface  of  the  propellant,  a  heat  balance,  as  suggested  by  the 
qualitative  description  of  the  combustion  zone  discussed  earlier,  indicates 
that  the  total  heat  received  by  the  surface  is  given  by 


-k  (jg)  =  F  +  Z1  P  exp  (kE^RI)  -  ffl  pr  (5) 
’x=0 

(1)  (2)  (3) 
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The  terms  on  the  right-hand  side  of  equation  5  represent: 

(1)  =>  heat  supplied  from  the  flame  zone  (gas -phase  reaction  zone  a 

finite  distance  away  from  the  surface) 

(2)  «  heat  generated  at  the  surface  by  heterogeneous  reactions 

between  the  solid  oxidizer  decomposition  gases  and  the  exposed 
fuel  matrix 

(3)  *  heat  released  because  of  gasification  of  propellant  at  the 

surface 


Before  the  oscillatory  response  to  induced  oscillations  can  be  studied 
with  the  preceding  equations,  the  effect  of  pressure  on  the  gas-phase  heat 
transfer  must  be  considered.  This  effect,  which  has  received  attention  in 
several  studies, ■**  »  is  the  result  of  a  complex  interaction 

between  diffusional  processes,  which  are  independent  of  pressure,  and 
chemical  reaction  processes,  which  are  highly  pressure  dependent. 


In  recent  studies,  Denison  and  Baum, ^  ^  Marxman, ^  ^  and  Friedly^ 

have  used  approximate  solutions  of  the  laminar  premixed  flame  equations  to 
characterize  the  gas  phase  contributions.  These  equations  generally  have 
the  form 


pr 


exp  (-Ef/RTf) 


(6) 


Unfortunately  the  use  of  equation  6  results  in  two  fundamental  limitations. 
First,  the  transient  analysis  based  on  this  equation  predicts  oscillations 
in  the  flame  temperature.  However,  thermodynamic  calculations  show  the 
flame  temperature  to  be  nearly  independent  of  the  pressure.  Under  the 
assumption  of  pseudoequilibrium  in  the  gas  phase,  the  flame  temperature 
should  also  remain  constant  under  oscillatory  conditions  as  well.  The  only 
reasonable  mechanism  for  a  shifting  flame  temperature  is  a  shifting  oxidizer- 
to-fuel  ratio,  but  such  provision  is  not  included  in  these  analyses. 

The  second  limitation  concerns  the  applicability  of  the  results  to  the 
experimental  data  obtained  in  the  present  studies.  Equation  6  does  not 
consider  the  heterogeneous  nature  or  "unmixedness"  of  the  composite  pro¬ 
pellant  combustion  zone.  Though  the  combustion  zone  of  double-base  propellants 
is  perhaps  one  dimensional,  this  certainly  is  not  the  case  for  composite 
propellants.  To  account  for  these  effects,  Summerfield''''  has  suggested 
that  the  characterizing  expression  shotild  be 


F"1  =  aP_1  +  bP~1/3  (7) 

(23) 

while  Fenner  has  developed  a  more  complex  expression  based  on  a  slightly 
different  approach  but  which  yields  the  same  approximate  result  as  equation 
7.  It  can  be  shown  that  equation  7  can  be  approximated  over  reasonable 
pressure  ranges  by 

1 

F  »  Fo  (P/P0)D  (8) 
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In  view  of  the  obvious  simplicity  of  equation  8  over  equation  7  and  the 
approximate  nature  of  equation  7,  equation  8  was  used  to  characterize  the 
gas  phase  effects  in  the  analytical  development. 

The  heat  generated  by  surface  reactions  has  been  incorporated  into  the 
surface  energy  balance,  equation  5,  by  assuming  the  overall  reaction  is 
first  order  in  pressure.  This  assumption  is  based  on  the  observations  of 
Anderson,  st  al,'-4'  which  showed  these  reactions  to  be  initially  controlled 
by  the  rate  of  absorption  of  reaction  species.  More  recently,  Hermance^° » ^ 
has  shown  theoretically  that  pressure-dependent  surface  reactions  can  con¬ 
tribute  significantly  to  the  overall  combustion  process  in  the  intermediate 
pressure  region  between  100  and  1,500  psi.  In  view  of  these  observations, 
it  is  reasonable  to  assume  that  the  effects  of  these  reactions  can  be 
incorporated  as  shown  in  equation  5. 

The  analysis  also  requires  an  expression  relating  the  regression  rate 
to  the  surface  temperature.  In  these  studies,  it  has  been  assumed  that  the 
linear  burning  rate  is  determined  by  some  overall  pyrolysis  rate  of  the 
solid  phase.  Mathematically,  this  can  be  expressed 

r  »  Z2  exp  (-E2/RTa)  (9) 

In  reality,  equation  9  represents  an  approximation  of  the  physically 
real  situation  in  which  several  reactions  are  occurring  simultaneously. 
Certainly,  a  more  elaborate  expression  could  be  used  where  the  burning  rate 
is  related  to  the  weighted  sum  of  a  series  of  Arrhenius  expressions.  How¬ 
ever,  the  hypergolic  ignition  studies  of  Anderson^'  have  shown  that 
significant  energy  release  can  result  from  interfacial  reactions  without 
substantial  gasification  of  reaction  products.  Thus,  to  a  first  approxi¬ 
mation,  one  would  expect  the  principal  effect  of  these  reactions  to  affect 
the  thermal  balance  at  the  surface  and  not  the  material  balance. 
Mathematically  therefore,  the  effects  of  these  reactions  have  been  included 
in  equation  5  but  are  not  explicitly  shewn  in  equation  9. 

To  complete  this  analysis,  the  temperature  deep  within  the  propellant 
remains  constant  at  the  initial  propellant  temperature.  Mathematically, 
this  can  be  expressed  as 


T  =  Too  at  x  ®  «  (10) 

Equations  4,  5,  8,  and  10  represent  an  approximate  one -dimensional  model 
of  the  combustion  zone.  The  equations  are  highly  nonlinear,  and  solutions 
require  finite  difference  methods  and  an  electronic  computer.  By  making 
certain  approximations  in  these  equations,  an  approximate  solution  for  the 
linear  response  characteristics,  i  e,,  the  acoustic  admittance,  can  be 
developed  to  guide  the  interpretation  of  experimental  data. 
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Baer,^  ^  in  his  analytical  studies  on  ignition  of  solid  propellants, 
has  shown  that  it  makes  little  difference  whether  the  energy  source  term 
is  in  the  boundary  condition  or  in  the  basic  differential  equations. 
Physically,  this  approximation  assumes  that  the  solid  phase  reaction  zone 
is  concentrated  near  the  propellant  surface  as  a  result  of  the  high  tempera¬ 
ture  gradient  in  the  solid  phase.  Because  the  temperature  gradient  is  high, 
this  observation  also  should  apply  to  analysis  of  combustion  instability. 
Thus,  as  a  further  approximation,  the  energy  source  term  in  equation  4  can 
be  eliminated,  which  yields 


(ID 


2.2.2  Linear  Response  Studies 
By  letting 


T(x,t)  «  T°  (x)  +  61  (x,t)  (12) 

where  0(x,t)  is  small  compared  to  T°(x),  equations  5,  8,  9,  and  11  become 

,  d6  dV 

7  r  =  —r  +  —  -  T°  0  (0)  (1-0)|1  exp  (-0  (13) 


ae^o.t) 


=  T°  (0).(q  n'  +T)  P/Pq  +  T°  (0).(rA-p)p 


(14) 


r0  (1  +  p) 

(15) 

(1  +  n'  P/P  ) 
o 

(16) 

Equation  13  can  the^be  solved  using  LaPlace  transforms  as  suggested  by 
Denison  and  Baum. '  '  When  the  boundary  conditions  given  by  equation  13 

are  applied,  the  result  becomes 


/GO  » 


_ dAsL±n-im _ 

(4j  +  4  (1-0) /s)  +r  A  0  +  4(l-0)/s 


(17) 
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The  acoustic  response  function  for  this  simplified  model  of  the 
combustion  zone  can  be  determined  by  investigating  the  oscillations  in  the 
burning  rate  produced  by  sinusoidal  pressure  oscillations  of  small  amplitude, 
that  is  to  pressures  having  the  mathematical  form 


P  =  P  (1+  fc) 

o  i 

(IB) 

Hence 

/(P)  =  e^Cs-a') 

(19) 

Substituting  equation  19  and  17  yields 

\k 

el 


(r+v)  “ 


|'(l-Wa^)  erfc  (-T  \Ja^  (lWa^)  erfc  (t 
- + - 


(Yx  -  aj_)  (Y2-  ai>  (Yx+  ax)(Y2+  ax) 


at 

e 


2Yi(Yi+1)  erfc  (-Y1‘/f)e  ^l"1^ 
+  — - - 

(Yx-  ax>  (Yj+  axXYx-Y^ 


2Yo(Y2+1)  erfc  (-Y?V?)e  (y2_1)t 

+  - - -  (20) 


(Y 2“  al>  <Y2+  ax)  (YX-Y2) 


where 

Y^  -  -2  (1+R)  ;  R  5 

yl  Y2  =  (1-2R  +  4S)  ;  S:^ 


(21) 

i 


The  steady  oscillatory  behavior,  provided  such  a  condition  exists  at  all, 
can  be  obtained  by  an  examination  of  equation  20  for  large  times. 


2c 2. 2.1  Stability  in  Self-Excited  Modes 

To  determine  if  steady  oscillatory  behavior  does  exist  for  a  particular 
set  of  combustion  parameters  (i.e.,  (P  -r  A  ) /\Jr  and  (l-0/\Jr) ,  the  behavior 
of  the  third  and  fourth  terms  in  equation  20  must  be  examined.  From  these 
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two  terms,  it  is  apparent  that  these  terms  will  disappear  with  time  to  produce 
a  steady  oscillatory  response  if  either 


Re  ( y  -  1)  <  0 


(22) 


or 

Re  (Y  ^  "  1)  >  0  but  erfc  (-y^  V7)-»  0  as  t  -*•  00  (23) 

From  the  properties  of  the  complimentary  error  function  having  complex 
arguments  and  equations  21,  22,  and  23,  the  regions  of  stable  oscillatory 
behavior  can  be  determined  as  a  function  of  the  parameters  (f3  -rA)/\jf  and 
(1-0) /ty.  The  resulting  stability  map  is  shown  in  figure  12  together  with 
the  relations  between  (f}-rA)/\}r  and  (1-0) /ty  separating  the  two  regions. 

Also  shown  are  the  regions  in  which  the  burning  rate  oscillations  increase 
without  limit,  either  in  an  oscillatory  manner  with  exponentially  increasing 
amplitude  or  exponentially  without  oscillation. 

This  figure  shows  that  stable  oscillatory  behavior  always  occurs  when 
p>f;,  ,  regardless  of  the  temperature  parameter  (1-0) /ty.  Hence,  if  the  net 
heat  of  pyrolysis  of  the  propellant  surface  (a  weighted  average  for  typical 
binders  and  ammonium  perchlorate  for  composite  propellants)  is  sufficiently 
endothermic  to  counter  the  effects  of  exothermic  pres sure- dependent  reactions, 
stable  oscillatory  behavior  will  occur.  If  not,  then  instability  in  the 
self -excited  modes  can  occur  (depending  on  (1-0) /\}r)  in  any  combustor.  All 
that  is  required  is  some  pressure  distribution  to  initiate  the  transient 
process. 

It  interesting  to  speculate  whether  in  fact  any  solid  propellants  have 
been  prepared  which  are  unstable  in  the  self-excited  mode.  It  seems  reason¬ 
able  to  assume  that  minor  pressure  disturbances  occur  in  any  combustor, 
including  the  normal  burning  rate  strand  bomb.  If  the  propellant  Is  unstable 
to  self-excited  modes,  these  disturbances  would  then  induce  burning  rate 
oscillations  of  increasing  amplitude  or  an  exponential  increase  in  burning 
rate  with  time.  These  burning  rate  disturbances  will  continue  until  second 
order  effects  somehow  limit  the  disturbances  or  until  the  propel! mt  Is 
consumed.  In  either  case,  it  seems  reasonable  to  expect  erratic  burning 
rate  data  at  the  very  best.  Thus,  the  principal  value  in  studying  self- 
excited  instability  is  that  it  provides  a  bound  on  the  combustion  parameter 
(1-0) /ij r  and  (p-TA)/t  which  cannot  be  violated  when  analyzing  experimental 
acoustic  response  function  data. 


2. 2. 2. 2  \coustic  Response  Function  ». 

Having  now  determined  the  regions  of  steady  oscillatory  behavior, 
equation  18  can  then  be  reduced  to  permit  prediction  of  the  acoustic  res¬ 
ponse  function.  To  perform  this  calculation,  equation  2  suggests  that  one 

should  use  a  complex  frequency  (i.e.,  CX  +  12  n  f  )  in  equation  20  to  obtain 

§  S 
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meaningful  comparison  between  theory  and  experiment.  This  would  account 
theoretically  for  the  effects  of  the  conditions  of  constantly  increasing 
amplitude  of  the  pressure  disturbance  which  are  used  experimentally  to 
evaluate  acoustic  response  functions.  However,  assuming  P  =  200  psi, 
c  ~  3,000  ft/sec,  p  =  0.06  lb/in.  ,  t  =  0.15  in. /sec,  and°Re  (p/e)  =  6, 
equation  2  indicates  ° 

a  /f  =  0.1 
S  S 

Therefore,  one  can  theoretically  study  the  response  to  fixed-amplitude 
pressure  oscillations  and  compare  the  results  to  experimental  data  obtained 
under  conditions  of  constantly  increasing  amplitude,  at  least  for  the  test 
conditions  and  propellants  studied  in  this  program.  Hence  to  =  i  27rfg  can  be 
substituted  into  equation  20  to  yield 


In  addition  to  predicting  the  steady  oscillatory  response  properties 
of  the  combustion  model,  equation  24  may  be  used  to  predict  the  theoretical 
pressure  exponent'  of  the  steady-state  burning  rate,  i.e.,  n  in  the  equation 

r  =  a  Pn  (25) 

o  o  v  ' 

From  equations  1  and  24 


n  «  Lim  Re  (p/e)  =  „  rA  j*  (26) 


These  results  indicate  the  acoustic  response  (p/e)  -  frequency  (X.)  relation¬ 
ship  is  a  function  of  three  parameters  .  \|f  4(1-0)  .  and  &-TA  or 

Kr+  q  n7'  \Jr  \Jf 

n,  4(1-0)  Tand/ft-r A  \  .  A  parametric  investigation  of  the  behavior  of 
t  \  ^  / 

equation  24,  for  the  stable  region  as  shown  in  figure  12,  is  shown  in  figure  13. 
A  number  of  interesting  observations  concerning  the  predicted  behavior  are 
apparent.  First,  the  calculations  indicate  that  high  ratios  of  maximum 
acoustic  admittance  to  burning  rate  pressure  expocant  occur  when  ({3-rA)/iji 
is  negative;  that  is  when  the  summation  of  all  surface  pyrolysis  reactions 
results  in  a  net  exothermic  process.  Under  these  conditions,  the  surface 
pyrolysis  reactions  which  occur  at  the  point  of  maximum  temperature  oscilla¬ 
tions  become  sources  of  energy  which  amplify  the  temperature  oscillations 
rather  than  sinks  which  attenuate  the  oscillations.  The  only  source  of 
thermal  damping  results  from  the  conduction  of  energy  into  the  solid  phase; 
as  a  result,  the  acoustic  response  function  is  greatly  increased. 
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Another  result  from  the  analysis  is  that  the  frequency  at  which  the 
maximum  admittance  occurs  depends  only  on  the  temperature  parameter,  S. 
This  is  apparent  both  by  differentiation  of  equation  20  with  respect  to  X 
and  from  the  parametric  calculations  shown  in  figure  13.  Hence  a  shift  in 
the  frequency  at  which  the  maximum  response  function  is  observed  can  only 
be  effected  by  a  change  in  S. 


It  should  be  noted  that  all  the  conditions  shown  in  figure  8  are  stable 
in  the  self-excited  modes.  Therefore,  a  propellant  can  be  stable  in  the  self- 
excited  modes  but  be  highly  unstable  in  the  combustion  chamber  as  a  result 
of  the  feedback  process  between  the  dynamic  pressure  response  of  the  combus¬ 
tion  chamber  to  burning  rate  oscillations. 


A  third  observation  concerns  the  effect  of  pressure- sensitive  surface 
reactions  on  the  predicted  behavior.  From  figure  13  it  is  apparent  that 
increasing  the  contribution  of  these  reactions  (i.e.,  increasing  r) 
significantly  increases  the  predicted  response  function.  This  then  raises 
the  question  as  to  the  influence  of  pressure  on  the  predicted  behavior. 

The  calculations  shown  in  figure  13  do  not  permit  the  investigation  of  the 
pressure  effect  directly  since  these  calculations  were  performed  by  holding 
the  parameters  ty,  {3,  and  0  constant,  which  in  effect  requires  that  rQ  and 
T°(0)  remain  constant  with  changing  P.  However,  changes  in  pressure  will 
obviously  change  both  rG  and  T°  (0).  To  circumvent  this  difficulty,  the 
dimensionless  groups  can  be  redefined  to  eliminate  the  surface  temperature 
as  follows: 


A  = 


w 

Z2CpE  2 


=  f  T  exp  (  ~~) 


(27) 


B  5 


F  R 
o 


Z2GpE2 


H  ir  e 


-1/* 


C?  R1  /E.  =  t0 
cn  2 


d  =  ah  r/c  e„  =  40 

v  p  2 


h'h 


Substituting  these  definitions  into  the  steady-state  energy  balance 
(equation  5  written  for  steady-state  conditions) 

B  exp  (1/t)  +  A  exp  (1-  A  )/4  =  +  D-C  (28) 
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Thereby,  the  value  for  the  dimensionless  surface  temperature  can  be  calcu¬ 
lated  for  each  specific  set  of  values  for  A,  B,  C,  D,  and  A  from  equation  28. 
Because  parameters  A  and  B  are  directly  related  to  pressure,  the  effect  of 
pressure  on  the  theoretical  acoustic  admittance  can  be  studied  using 
equations  24,  27,  and  28.  The  results  are  shewn  in  figure  14  for  one  partic¬ 
ular  case  where  surface  reactions  contribute  approximately  20%  of  the  total 
energy  flux  to  tlje  solid  surface.  It  should  be  noted  that  the  dimensionless 
frequency  4ncu/r  contains  the  burning  rate  at  the  pressure  of  interest  and 
hence  contains  a  pressure  effect  which  is  not  apparent  in  figure  14  .  To 
remove  the  effect,  the  frequency  has  been  renormalized  with  respect  to 
instead  of  f .  Thus,  the  modified  dimensionless  frequency  becomes 

X  '  =  ka  co/Z*,  (29) 

and  is  directly  proportional  to  the  dimensional  frequency w.  On  this  basis, 
the  results  shown  in  figure  14  have  been  replotted  are  are  shown  in  figures 
15  and  16. 

From  figures  14  through  16,  it  is  apparent  that  the  model  predicts  the 
frequency  at  which  the  maximum  acoustic  response  increases  with  increasing 
pressure.  This  results  from  the  pressure  effect  on  the  burning  rate  and 
the  resultant  effect  on  the  characteristic  time  of  heat  conduction.  This 
effect  should  be  universal  to  all  combustion  models.  These  figures  also 
show  the  maximum  admittance  and  pressure  exponent  increase  as  the  pressure 
increases.  The  magnitude  of  all  these  effects  depend  on  the  values  of  the 
various  combustion  parameters,  and  further  parameter  studies  are  required 
to  determine  which  parameters  are  the  most  significant.  These  studies  are 
currently  in  progress. 


2.2.3  Comparison  of  Experiment  and  Theory  in  the  Linear  Region 

When  the  theoretical  calculations  are  considered  in  the  light  of 
experimental  observations,  some  interesting  correspondence  and  some  disturbing 
differences  can  be  noted.  First,  if  the  analysis  does  approximately  model 
the  actual  combustion  process,  then  the  high  ratios  of  the  maximum  acoustic 
response  function  to  the  pressure  exponent  of  the  burning  rate  signify  the 
net  surface  pyrolysis  processes  are  either  thermally  neutral  or  are  exothermic, 
i.e.,  R  is  zero  or  negative.  Unless  this  is  the  case,  the  model  cannot 
account  for  the  Re  (p/e)max/n  values  of  4  to  7  which  are  observed  experi¬ 
mentally.  This  would  suggest  the  exothermic  decomposition  of  ammonium 
perchlorate  plays  an  important  role.  This  observation  is  at  least  consistent 
with  the  theoretical  calculations  of  Hermance^  '  and  with  the  experimental 
data  obtained  by  Waesche  et  al.',/°' 

If  R  is  negative,  the  shift  in  frequency  at  which  the  maximum  response 
occurs  when  coated  oxidizer  is  used  suggests  the  coatings  reduce  S  by 
reducing  the  overall  activation  energy  of  the  surface  pyrolysis  reactions 
or  increase  T°(0).  A  reduction  in  reduces  both  the  parameters  R  and  S. 
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Figure  14.  Acoustic  Response  Function  —  Moderate  Surface  Reaction 

Term  Effect  on  Pressure 
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Figure  16.  Predicted  Acoustic  Response  for  n'=1.0  as  a  Function  of 
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IX  R  is  then  negative,  figure  13  indicates  that  both  the  maximum  response 
function  and  the  frequency  at  which  this  maximum  occurs  are  reduced.  This 
behavior  was  observed  experimentally  with  the  Kel-F  and  Viton  A  coated 
oxidizer. 

However,  the  effect  of  other  coatings  was  to  shift  the  maximum  response 
frequency  but  to  not  change  the  maximum  response  itself.  In  the  context  of 
the  analysis,  such  behavior  can  only  result  from  a  compound  effect  such  as 
a  concurrent  decrease  in  E  and  an  increase  in  either  q,  n',  or  T.  The 
changes  in  the  chemistry  which  effect  this  change  are  not  apparent,  however. 
Clearly  changing  -he  activation  energy  suggests  a  shift  in  the  contribution 
of  the  various  reactions.  One  might  then  expect  a  concurrent  shift  in 
burning  rate  and  surface  temperature,  but  these  secondary  changes  are  not 
apparent  from  the  analysis.  More  analysis  of  the  results  shown  in  figure  13 
are  required  to  elucidate  the  changes  produced  by  the  coatings  within  the 
context  of  the  model. 

The  predicted  effect  of  pressure  shown  in  figure  14  is  consistent  with 
experimental  observations  of  Ibiricu'1  '  and  suggest  that  pressure-dependent 
surface  reactions  cannot  be  a  major  contribution  to  the  overall  combustion 
process,  at  least  at  low  pressures  for  this  propellant.  Under  conditions 
of  even  moderate  pressure-dependent  reactions  at  the  surface,  the  predicted 
burning  rate  pressure  exponent  increases  with  increasing  pressure.  The 
predicted  effect  contrasts  experimental  burning  rate  data  for  the  propellants 
studied  and  thereby  raises  the  question  of  the  validity  of  the  basic  model. 

It  should  be  noted  that  pressure  exponents  are  generally  observed  to  increase 
significantly  at  elevated  pressure,  i.e.,  2000  to  5000  psi  for  many  composite 
propellants. 

The  predicted  effect  of  pressure  on  the  acoustic  response  function 
shown  in  figures  13  through  16  has  some  similarity  with  the  data  reported 
in  figures  5  and  6  and  with  the  experimental  observations  of  Ibiricu. 

Both  the  magnitude  of  the  maximum  acoustic  response  functions  and  the  fre¬ 
quency  at  which  the  maximum  occurs  increase  with  increasing  pressure  at  lower 
pressures.  However,  there  are  some  serious  differences  between  theory  and 
experiment.  At  higher  pressures  the  maximum  response  function  continues  to 
increase  with  pressure  while  Ibiricu's  data  show  the  reverse  effect.  Further¬ 
more,  the  predicted  shift  in  the  frequency  of  the  maximum  response  is  much 
larger  than  observed  experimentally.  Also  the  predicted  magnitude  of  the 
maximum  response  functions  are  substantially  less  than  the  experimental 
data.  This  indicates  that  the  combustion  parameters  chosen  for  the  para¬ 
metric  study  are  not  characteristics  of  the  propellants  studied.  The 
studies  currently  in  progress  are  directed  towards  determining  the  approximate 
combustion  parameters  for  these  various  propellants  to  establish  if  realistic 
comparisons  between  theory  and  experiment  are  possible. 


2.2.4  Nonlinear  Response.  Studies 

Nonlinear  analytical  studies  have  been  directed  toward  two  aspects  of 
the  acoustical  behavior:  (1)  the  effect  of  subsurface  reactions  and  (2) 
the  effect  of  second  order  effects  in  the  absence  of  subsurface  reactions. 
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In  the  first  technical  summary  report,^  ^  the  equations  describing  the 
effects  of  subsurface  reactions  were  developed  together  with  the  finite 
difference  algorithm  for  their  numerical  solution  on  the  digital  computer. 
This  development  is  also  presented  in  appendix  B. 

Both  the  computer  code  for  predicting  the  steady-state  temperature  dis¬ 
tribution  within  the  solid  and  the  code  to  predict  transient  behavior  have 
been  developed  and  successfully  tested.  Good  numerical  stability  and 
convergence  were  demonstrated. 

However,  difficulty  arises  when  the  results  from  the  steady-state  code 
are  incorporated  into  the  transient  analysis.  Unfortunately,  the  calculated 
profile  is  never  precisely  equal  to  the  exact  steady-state  profile  because 
infinite  computer  time  would  be  required.  The  transient  analysis  considers 
these  errors  as  heat  fluxes  from  some  unknown  source  or  sources.  These 
numerically-generated  fluxes  are,  of  course,  erroneous  but  unfortunately 
they  are  approximately  of  equal  magnitude  with  the  transient  fluxes  generated 
by  the  oscillating  pressure.  The  calculations  conducted  to  date  have  shown 
that  these  errors  cause  severe  convergence  problems.  The  only  apparent 
method  to  circumvent  this  problem  is  to  use  similar  mesh  spacings  in  both 
the  space  and  time  variables.  However,  because  the  computer  costs  become 
prohibitive  when  the  required  mesh  spacings  are  used,  further  studies  along 
this  avenue  were  terminated. 

The  principal  justification  for  using  numerical  techniques  was  that 
subsurface  reaction  effects  on  both  the  linear  and  nonlinear  behavior  could 
be  determined.  If  these  effects  are  assumed  to  be  small,  then  analytical 
techniques  can  be  used  as  an  alternate  method  to  estimate  the  nonlinear 
behavior  of  the  combustion  model.  To  apply  this  approach,  temperatures 
within  the  solid  and  pressure  were  expanded  in  the  form 


p  =  P  fl  +  e'  e  exp  (Da)  +  e'  e  exp  (Da)l  (30) 
°  l  o  i 


r(t)  =  r 


T  (x,t)  = 

T(x,0)  +  e'  0^  (x,t)  +  e' 

e2  (x,t) 

r 

,2  2 

7 

1  +  e,01(O,t)E2  +  e  /E201(°,t)\  + 

£'  e2(o,t)E2 

RT2  (0,0)  2  \rT2(0,0)  / 

rt2(o,o) 

(31) 


(32) 


Equations  30,  31,  and  32  were  then 
through  11  and  collecting  terms  of 


substituted  into  equations  5  and  8 
the  order  ^  yields 
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ff2  =  *X  +  +j 

'^(O.t)’ 

2 

dT(x.O)  ei(0,t)  .  9Ci(0>t) 

3t  d£2  d£  i 

*T(0,0)  2 

.n(  o,o). 

d£  *T(0,0)  9^> 

(33) 


at  E,  =  0 


1.^ 


d0o(x,t)/T(O,O)  =  Oin+F)  e,(t)  +  (F  A-ft)  +  (T  A  -ft) 


4T(0,0) 


0^0,  t) 


n(o,o) 


I  nn'(n'-l)  2,,,  ,  rA  ¥°- ■*  1< 

2  €  1  <  4, 1(0,0) 


As  £  -*  °o  @2  -*■  0 


(34) 

(35) 


(14) 

Friedly  has  noted  that  the  real  components  in  the  nonhomogeneous 
terms,  such  as  0^  (0,0)  should  be  products  of  the  real  components.  In  other 
words,  these  cross  products  must  be  regarded  at  Re  (0^)  •  Re  (0^)  when  the 
complex  notation  is  considered.  This  observation  then  leads  to  the  conclu¬ 
sion  that  these  nonhomogeneous  terms  have  the  form 


Re 


(xx)  •  Re  (x2)  =  x12  (£)  +  x12  (Q  exp  (iXx/2) 


(36) 


and  hence  ©2  must  have  the  form 


62  (0  »  (0  +  ec2  (0  exp  (i.\x/2) 


(37) 


Thus,  consideration  of  second  order  effects  leads  to  two  solutions,  one 
which  is  independent  of  time  and  one  which  oscillates  at  the  second  harmonic 
frequency.  The  time -independent  solution  represents  a  shift  from  the  steady 
state  and,  therefore,  is  the  more  important  second  order  effect. 

Substituting  equation  36  and  37  into  equations  33  and  34  and  solving 
for  the  time-independent  component  yields 


d2p  °  (D  d  ,x°  (0 

- -4 - +  - - -  -S 


(0)  + 


Re2(u/e)+Im2(u/e) 


+  Re  WO  ±  l.ni  illM  .  E1  cos(MiQ  -  p1  sin^))  exp(Mr0  =  0  (38) 
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Substituting  equation  41  and  24  into  equation  32  therefore  permits  the 
change  in  burning  rate,  i.e., 


The  results  of  parameteric  studies  conducted  to  date  are  shown  in  figures 
17  and  18  for  the  same  cases  shown  in  figures  14  through  16.  In  figure  17, 
the  effect  of  increasingly  exothermic  pyrolysis  reactions  is  seen  to  increase 
the  magnitude  of  the  shift  in  steady- state  burning  rate.  Comparing  the 
results  on  figure  14  with  those  on  figure  17  shows  the  maximum  acoustic 
response  and  the  maximum  decrease  in  burning  rate  both  occur  at  the  same 
frequency. 

It  should  be  noted  that  the  maximum  change  shown  in  figure  17  still 
represents  a  small  shift  in  burning  rate.  Typical  values  for  e^,  in  the 
T-burner  are  0.1  to  0.3  and  hence  the  maximum  decrease  in  burning  rate 
shown  in  figure  16  is  approximately  TU.  This  value  is  of  the  same  order  of 
magnitude  as  the  experimental  observations  of  Eisel.'-  '  Since  there  is  no 
assurance  that  the  constant  parametric  values  used  to  prepare  Figure  16  are 
characteristic  of  Eisel's  propellants,  the  differences  between  theory  and 
experiment  have  no  significance. 

Figure  18  shows  the  predicted  effect  of  pressure  on  the  burning  rate 
changes,  again  for  the  same  basic  set  of  parameters  used  in  figures  14 
through  16.  Two  results  are  of  particular  interest.  First,  there  appears 
to  be  little  effect  of  pressure  on  the  shift  in  burning  rate.  This  holds 
for  both  the  case  where  n'=1.0  as  shown  and  for  the  case  of  n'=0.5.  In 
view  of  the  larger  pressure  affect  on  the  maximum  acoustic  response  function 
when  n'=0.5,  it  appears  that  the  second  order  pressure  effect  is  generally 
small.  However,  this  should  be  taken  as  a  preliminary  observation  since  a 
wide  range  of  combustion  parameters  have  not  been  considered. 

The  second  result  shown  in  figure  18  is  the  very  small  increase  in 
burning  rate  at  the  lower  frequencies  under  oscillating  pressure  conditions. 
This  result  is  real  when  the  limits  of  equation  44,  41,  and  24  are  considered 
and  results  from  the  particular  set  of  combustion  parameters  chosen  for 
study.  This  effect  appears  to  be  extremely  small  in  any  event  and  is  there¬ 
fore  more  of  a  theoretical  curiosity  than  of  a  practical  interest. 


Figure  17.  Effect  of  D  on  Predicted  Shift  in  Burning  Rate 


Figure  18.  Effect  of  Pressure  on  Predicted  Shift  in  Burning  Rate 
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APPENDIX  A 

EXPERIMENTAL  DATA 


Data  for  the  propellants  'tested  on  the  present  program  are  tabulated 
in  this  appendix.  In  addition,  the  acoustic  response  data  are  shown  as  a 
function  of  the  measured  growth  frequency.  Previously  reported  results^^) 
were  made  using  a  frequency  based  on  8  it  f  /r2  .  This  latter  form  was 
suggested  by  the  theoretical  studies. 
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Propellant:  /C/~y  ?3Z>/ 

Propellant  density:  <7£>6>3  /£*,/£»+  ^ 
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Propellant:  CC7/7  ?S?>/ 

Propellant  density:  ^ 

Velocity  of  sound  in  product  gas  (c): 
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Propellant:  ^XX  -  3£2>  Z  /?/?  CL&0.  OV  as  ^ /s'%  A'<<X'/r'  e 

Propellant  density:  0 '  OX'S  ^  °  &  ^"3 

Velocity  of  sound  in  product  gas  (c) :  3 4,000 
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Propellant  density:  <2-0  63  £ 
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Propellant:  U7Y-  ?S  76  7  ( 6?fl7 bef  cos-//  /%  <9/47?-0- 
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Velocity  of  sound  in  product  gas  (c)  :  3>  ^> ,  5 --Wv,  /^rC-c.  . 
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Propellant: 

Propellant  density:  &.  &S?  ^ 

Velocity  of  sound  in  product  gas  (c)  :  £2>0  /k^^r 


FIGURE  I.  EFFECT  OF  COATING  OR  ADDITION  TO  PROPELLANT  MIX 


FIGURE  2.  EFFECT  QF  KEL-F  COATINGS 


FIGURE  3.  EFFECT  OF  VARIOUS  COATINGS 


FIGURE  4.  EFFECT  OF  KEL-F  AND  VITON  COATINGS  ON  ACOUSTICAL  RESPONSE 


FIGURE  6.  EFFECT  OF  AP  PARTICLE  SIZE 


FIGURE  7.  UTREZ  BINDER  FOR  MOVIES 


UTC  2136-TSR2 


71 


UTC  2136-TSR2 


Exact  solution  of  equations  18,  27,  and  28  is  difficult  because  of  the 
nonlinearity  of  these  expressions.  In  order  to  circumvei."  this  problem,  a 
finite  difference  algorithmn  was  developed  for  digital  com  ^‘•er  solution. 
Equation  B-5  in  continum  form  replaced  the  finite  differe’  >  ipproximation 
as  shown  by  equation  L-7  where  the  mesh  points  of  the  difi  ce  scheme  are 
represented  by  i  from  i  =  1  at  the  propellant  surface  to  a  >  point  nAt, 
within  the  solid  propellant. 


Gas  Phase 


Solid  Phase 


'+  +  +  + 
2  3  4  5 


For  i  =  2,  3,  4,  i,  n 


d0_±  I B  +B.e-E/RTj? 

3t  /  +  ro  2AC  +  B  6  1  ^  ' 

Writing  the  boundary  condition  equation  at  £  =  0  in  finitie  difference 
form  results  in 


0-40+30 
8  0  _3 _ 2 _ i 

ar,  “  2M 


1  -E  /RT^  0  a  .  .  . 

-m 


(B-8) 


This  equation  may  be  differentiated  with  respect  to  the  time  variable,  r , 
to  give 


9  0  4  90  9  0  \  /  .  -  -E  /RT  0 

— i _ -  +3  — -  )  /  2  A  4  =  ill/^  +  M  e  00  1 

9t  9t  9t/  /  9r  9r 


0“  RT 
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.  2  00 


OW  n  /  C. 
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(B-9) 
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Rearranging 
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(B-10) 


Thus,  equation  B-7  plus  the  n-1  equations  of  the  form  of  equation  B-4 
provides  a  set  of  n  simultaneous  differential  equations  which  can  be  solved 
with  the  aid  of  a  digital  computer.  From  equations  8  and  18  we  obtained 
P,  3P/3r  and  3v|/3t  and  the  time  differential  form  of  equation  9. 

In  these  studies  a  method  of  solution  has  been  used  which  is  quite 
unique  in  solving  the  transient  diffusion  and  chemical  reaction  equations. 
First,  a  system  of  ordinary  differential  equations  is  generated  from  the 
relevant  partial  differential  equation  by  discretization  in  terms  of  the 
space  variable,  and  then  the  sequence  of  values  of  the  dependent  variable  at 
the  mesh  points  are  regarded  as  new  dependent  variables.  This  leads  to  a 
system  of  ordinary  differential  equations  (in  derivatives  with  respect  to 
time).  This  means  of  solution  of  boundary  value  problems  was  extensively 
exploited  for  some  tiim.  in  analogue  computer  work  for  solution  of  such 
problems  as  transient  conduction  of  heat  in  solids.  In  analogue  computer 
work,  the  limitation  is  usually  the  capacity  of  the  computer.  The  relatively 
recent  development  of  digital  computer  subroutines  ("package  solvers") 
"analogue  discretization"  methods  to  a  competitive  means  for  solving  certain 
types  of  partial  differential  equations. 

Tne  specific  advantages  of  this  approach  to  solution  of  boundary 
problems  may  be  listed  as  follows: 

A.  The  programming  time  is  reduced  to  a  minimum  as  most  of  the 
difficult  programming  has  been  done  in  the  development  of 
the  ordinary  differential  equations  subroutine.  Also,  the 
computer  running  time  for  the  differential  equations  approach 
seems  to  offer  an  advantage  compared  to  problems  involving 
differencing  in  both  the  space  and  time  steps. 

B.  Numerical  convergence  problems,  i.e.,  convergence  problems 
arising  in  numerical  solution  of  the  partical  differential 
equations,  are  relegated  to  problems  of  convergence  of  solu¬ 
tions  of  systems  of  ordinary  differential  equations.  Most 
subroutines- for  solution  of  ordinary  differential  equations 
have  some  means  of  assurance  of  numerical  convergence  of  the 
solution . 
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C.  Conversion  of  a  boundary  value  problem  into  a  system  of 
ordinary  differential  equations  provides  a  means  of  handling 
nonlinear  difficulties.  This  point  is  particularly  important 
because  many  of  our  equations  involve  an  exponential  Arrhenius 
term  as  well  as  reactant  product  terms. 

D.  The  capacity  of  the  differential  equation  "package  solvers" 
is  sufficiently  large  to  allow  use  of  a  fairly  large  number 
of  mesh  points  in  the  discretization. 

The  study  of  the  acoustic  response  requires  that  the  oscillations  occur 
about  the  steady  state.  Hence,  the  analysis  must  first  predict  the  steady- 
state  conditions  for  a  given  set  of  combustion  parameters  and  then  calculate 
the  transient  behavior.  To  do  this,  two  approaches  are  possible.  One 
approach  is  to  assume  some  initial  temperature  distribution  within  the  solid 
and  let  the  transient  equations  relax  to  the  appropriate  steady  state.  The 
pressure  oscillations  are  then  started  because  the  steady  state  is  reacted. 

The  second  approach  is  to  neglect  the  transient  terms  and  calculate 
the  steady-state  conditions  by  r.  trial  and  error  process.  In  this  process, 
the  surface  temperature  is  assumed,  and  the  equations  are  solved  using  a 
forward  differencing  procedure  for  the  temperature  deep  within  the  solid. 

The  surface  temperature  is  then  adjusted  until  the  desired  temperature  in 
the  propellant  is  obtained.  The  resulting  steady-state  profile  is  then  used 
in  the  transient  model  and  the  oscillations  initiated  immediately. 

Of  these  two  possible  approaches,  the  second  is  preferred  for  economic 
reasons.  The  calculations  using  the  forward  difference  approach  can  be  done 
rapidly  on  the  computer  in  comparison  to  the  transient  calculations.  Hence, 
this  is  the  approach  used  in  these  studies. 

As  mentioned  in  the  main  body  of  *:his  report,  even  this  approach  did 
not  prove  to  be  economical.  The  difficulty  arose  when  the  steady-state 
profile  was  inserted  into  the  transient  equations.  Even  though  close  con¬ 
vergence  was  specified,  the  predicted  steady  state  never  matches  the  exact 
steady  state.  When  such  information  is  inserted  in  the  transient  equations, 
further  relaxation  to  the  true  steady  state  occurs.  This  relaxation  occurs 
during  the  initiation  of  the  oscillations  and  this  compounding  effect  was 
found  to  require  extremely  small  ratios  of  At /At,2  to  maintain  numerical 
stability.  This,  in  turn,  results  in  extremely  long  computer  running  time 
(i.e.,  1/2  to  1  hr)  before  the  steady  oscillatory  behavior  is  achieved. 

Many  such  calculations  would  be  required  to  acess  the  effects  desired, 
therefore,  these  studies  were  terminated  because  of  the  excessive  cost. 
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